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Electric response of spin polarization in two-dimensional electron gas with structural inversion 
asymmetry subjected to a magnetic field was studied by means of the linear and non-linear theory 
and numerical simulation with the disorder effect. It was found by Kubo linear reponse theory that 
an electric resonant response of spin polarization occurs when the Fermi surface is located near 
the crossing of two Landau levels, which is induced from the competition between the spin-orbit 
coupling and Zeeman splitting. The scaling behavior was investigated with a simplified two-level 
model by non-linear method, and the resonant peak value is reciprocally proportional to the electric 
field at low temperatures and to temperature for finite electric fields. Finally numerical simulation 
illustrated that impurity potential opens an enegy gap near the resonant point and suppresses the 
effect gradually with the increasing strength of disorder. This resonant effect may provide an efficient 
way to control spin polarization by an external electric field. 

PACS numbers: 



I. INTRODUCTION 

Recently electric or non-magnetic generation of spin 
polarization in semiconductors has attracted a lot of 
interests because of its potential application in spin- 
tronic devices of semiconductors^. In this attempt 
the motion of electron interacts with its spin via the 
spin-orbit coupling, which provides a possible way to 
control electron spin by electric fields. Several exper- 
iments were resolved to produce spin polarization in 
quantum wells^££, strained semiconductors^, and hole- 
doped hetero-j unction 8 . In a two-dimensional electron 
gas (2DEG) with structural inversion asymmetry, it was 
understood that the spin-orbit coupling generates an ef- 
fective momentum-dependent field to induce a net bulk 
spin polarization by electric fields or current o 9 ' 10 ! 11 ' 12 . 
Another mechanism has also been proposed to gener- 
ate spin polarization in a bulk 2DEG in the presence 
of in-plane magnetic fields and electric fields^. Electric- 
induced spin accumulation near the boundary of sam- 
ple was already observed experimentally in either n- or 
p-doped semiconductors as a consequence of spin Hall 
effec t 14 ' 15 i 16 . The ac-field method has also been applied 
to induce and detect spin polarization. The spin orienta- 
tion was achieved by the excitation of a high-frequency 
electric fielcUi. An ac-field detection of spin resonance 
was also discussed extensivel y 18 i 19 . 

Generally speaking, we may introduce the electric spin 
susceptibility \e t° describe the response of spin polar- 
ization S a to an external electric field Ep2£ 

S a {E) = S a {Q) + xfE p , (1) 

where ^(O) is the spin polarization in the absence of 
electric field. Usually the electric spin susceptibility is a 
tensor, not a vector. The spin polarization is determined 
by the symmetry of spin-orbit coupling of the system. 
To generate spin polarization efficiently, a large electric 
spin susceptibility is expected. In this paper we pro- 
pose an electric resonant response of spin polarization in 



2DEG with Rashba spin-orbit coupling. First the linear 
response theory shows that the electric spin susceptibility 
becomes divergent when the crossing point of two Lan- 
dau levels is accidentally located near the Fermi surface. 
The additional degeneracy of two Landau levels is at- 
tributed to competition between the spin-orbit coupling 
and the Zeeman splitting. Then a simplified two-level 
model was proposed to investigate the scaling behaviors 
of the resonant electric spin susceptibility. The resonant 
values of the electric spin susceptibility decay with either 
the applied electric field at low temperatures or with the 
temperatures in a weak electric field. Finally we take 
into account the disorder effect, and apply the truncation 
approximation to study the syetem numerically, which 
is beyond the linear response theory. The dependence 
of the divergent behavior on the electric field and tem- 
perature was presented for finite disorder strengths. It 
is damped with the electric field and temperature when 
the energy scale of electric field and temperature is larger 
than the impurity potential. The numerical consequence 
is consistent with the result of linear response theory. 



II. MODEL HAMILTONIAN AND LINEAR 
RESPONSE THEORY 

We start with a 2DEG with Rashba coupling confining 
in a two-dimensional plane L x x L y . The model Hamil- 
tonian in the presence of a perpendicular magnetic field 
B is given by 

H = H Q + H R , (2) 

where 

#o = TT—^n 2 - -g s Ms B a z , (3) 
2m* 2 

and the Rashba coupling 
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Hr = -r(Jl x CTy - U y a x ). (4) 

m* is the effective mass of the electron, II = P+|A is 
the kinetic momentum, g s is the Lande-<? factor, and [Ib 
is the Bohr magneton. We take the Landau gauge A = 
y Bx and consider the periodic boundary condition in 
x-direction. The Hamiltonian can be solved analytically 
with the eigenvalues 



e„, s = hw c \n + -^/(l - g) 2 +8nr] R J (5) 

where n is a non- negative integer and s = dz l 21 ' 2 , 
Here the cyclotron frequency lo c = eB/m*c, the mag- 
netic length lb = \J hc/eB, the dimensionless g-factor 
g = g s m* /2m e and effective coupling tjr = a/lbhuj c . The 
corresponding eigenvectors are expressed as 

\nks) = ( lc ^ nk \ , (6) 

V - sc ns<Pn-lk J 

where </) n k is the eigenstate of the n-th Landau level with 
k = p x the good quantum number because of the pe- 
riodic boundary condition along the x-direction. Each 
Landau level has a degeneracy = L x L y /2irl 2 . c^ s are 
the coupling parameters of spin up and down. c^ +1 = 1, 

c^_! = 0, and c± s = 1/ \J 1 + (u n =F Sy/1 + u 2 ) 2 with 

u n = (1 — g)/V&nrjR for n > 1—. One of the remark- 
able features of the spectra is the additional crossing of 
Landau levels, which is generated by the competition 
between Rashba coupling and Zeeman splitting such as 
£n.s=+i = £n+i,«=-i if the integer n and the magnetic 
field B satisfy 

^(1 - g) 2 + 8 n r? R +y/(l 9? + 8 (n + 1) ^ = 2. (7) 

One important factor in this system is the filling factor, 
i.e., the ratio of number of charge carriers to the Lan- 
dau degenercay v = N e /Nk = 2i:l 2 n e (n e is the density 
of charge carriers). For a specific density, the filling fac- 
tor is proportional to 1/B, and for a specific field, it is 
proportional to the density n e . 

Now we come to study the electric response of spin 
polarization when the Fermi surface is located near the 
resonant point. We apply an in-plane weak electric field, 
say, E y along y-direction, then the electric spin suscepti- 
bility can be evaluated by means of the Kubo formula in 
the weak field limit 2 ^, 

^ay _ e ^ -j-^ (fn's' ~ fns) 

" nn'kk'ss' ' 

(Kubo) 

x (nks\ S a \n'k's') (n'k's'\ v y \nks) 

where v y is the velocity in y-direction, f ns are the Fermi- 
Dirac distribution function and r is the lifetime of the 
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FIG. 1: Xe versus & t the temperature T — for 
weak electric fields by means of the Kubo formula ?? 
when t — ► oo. The parameters are taken as: 
a = 4.16 x Kr n eVm, 
n e = 0.118 x KT 16 m- 2 ,m* = 0.05 m e , g s = 4. The unit 
of x v E v is h/inl 2 b N/C. 



system. The denominator in the formula indicates that 
the electric spin susceptibility may become singular when 
two energy levels are degenerate near the Fermi surface 
and for a long lifetime r, i.e., Xe ma y become divergent. 
For a finite r or temperature it is still finite. 

For the purpose of numerical calculations, without loss- 
ing of generality, we take the model parameters such that 
the lowest energy level |1) = \n = 0, k, s = +1) and the 
first excited level |2) = \n = 1, k, s = —1) cross at a crit- 
ical magnetic field B = 9.8T. We plot Xe versus the 
inverse of the magnetic field 1/B in Fig. [T]for a fixed elec- 
tron density n e . It is found that a resonant peak occurs 
at the critical magnetic field corresponding to a filling 
factor v = 0.5. The peak indicates that it is possible for 
a weak electric field to induce a finite spin polarization 
at the resonant point. 

It is worth stressing that the Kubo formula is a result of 
perturbation and is valid only when E — > (but E 0). 
At a first glance the divergence of spin susceptibility in 
Fig. [T] may be unphysical because it might be caused by 
the approaximation of perturbation in the linear response 
theory. To clarify the problem, we have to go beyond the 
linear response theory and investigate the problem by 
means of the non-perturbative approach. 



III. NON-LINEAR BEHAVIORS IN SIMPLIFIED 
TWO-LEVEL MODEL 

In order to understand the physical origin of the reso- 
nance and to analyze the scaling behaviors of Xe a * or 
near the resonant point, we consider a simplified two-level 
model around the resonant point. This is based on the 
fact that the dominant contribution to the electric-field- 
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induced resonance is attributed to the energy crossing 
of the two levels. This truncation approaximation re- 
duces an infinite-dimension problem into a two-level one. 
As the two-level problem can be solved analytically, the 
effect of electric field can be taken into account. The 
method goes beyond the perturbation method and the 
nonlinear behaviors of spin polarization to external field 
can be revealed. For this purpose we are only concerned 
with the two nearly degenerated levels |1) and |2) and 
ignore the contribution from other levels of higher ener- 
gies. Then the total Hamiltonian including the electric 
potential of external field can be reduced to a 2 x 2 ma- 
trix: 



H* = 



S VE 

v E -e 



(8) 



where e = \ (ei,_i - £o,+i); and v E = (l\eE y y\2) = 
r)R eElbCQ +1 cf _i The off-diagonal element induced by 
the electric potential lifts the degeneracy and opens an 
energy gap 2\ve\ at the resonant point of e = 0. In Fig. 
[H the energies of the two levels and the spin polariza- 
tions for the lower energy states are plotted with repect 
to the inverse of magnetic field. In the absence of electric 
field the two levels cross at the resonant point B = 9.8T 
and each state has almost opposite spin polarization in 
z-direction. From the point of view of the lowest en- 
ergy state, the spin polarization S z has a jump near the 
crossing point as shown in Fig. [2] (b). In the presence of 
electric field, the two states of |1) and |2) will be admixed 
due to the off-diagonal term in H * . The term will lift the 
degeneracy of the two levels near the resonant point. In 
the lowest energy state, the spin polarization rotates in 
the z-y plane from positive spin S z to the negative S z 
when the magnetic field sweeps over the resonant point. 
At the resonant point, the system is almost polarized 
along the y-direction in this two-level problem. Compar- 
ing with the case in the absence of electric field, we notice 
that non-zero S y is induced by a weak electric field. Be- 
cause there is a finite density of states at this point due 
to the Landau degeneracy, the total spin polarization in 
y-direction becomes finite. As a result the electric spin 
susceptibility Xe becomes divergent at the point. 

If we denote by the expectation value of spin po- 
larization S y of the two split levels, the overall response 
of the spin polarization to a finite external field can be 
calculated by: 



Xe = 



i (Sy"U + s;.f-) E -(s+u + s-.f-) E=0 



E 



(9) 

with f± the Fermi-Dirac distribution of the two levels 
and /_ + /+ = v. With other parameters the same as in 
Fig. ? , electric spin susceptibility Xe? near the resonant 
point is plotted as a function of 1/B in Fig. [3] (a) for 
different electric fields. From Eq. @ we can rewrite it 
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FIG. 2: (color online) (a) and (b): Energy spectra for 
the two degenerate levels in the absence of the electric 

field and the correponding spin polarization of the 
lower-energy states, (c) and (d): Energies for the two 
anti-crossing levels in the presence of the electric field 
and the corresponding spin polarization for the 
lower-energy level. The higher energy states are denoted 
by dotted line and the lower are denoted by solid line. 
The arrows indicate the orientation of S z . The spin 
polarization S z and S y for the lower-energy states are 
denoted by blue and red line respectively. The energy is 
in unit of hoj r . The parameters are taken the same as in 



Fig. 1: a = 4.16 x 10" n eVm, 



0.05 m e , g s = 4. 



in the form: 



yy _ /-(I - /-+ 
Xe TTl^E 



VE 



VE 



> + \var/kT 



(10) 

Near the resonant point e — > 0, we have Xe ^ V-^ a ^ 
low temperatures kT -C \ve \ , meanwhile Xe oc 1/T for 
a weak electric field kT 3> |«e| . It depends on the energy 
scales of the electric field energy and temperatures. In 
Fig. [3] (b) and (c) the dependence of the peak value of 
x'e on the electric field and temperature are plotted. A 
finite spin polarization can be induced by a weak electric 
field, but x'e decays with the electric field and the tem- 
perature. This result is consistent with the result of the 
Kubo formula in Fig. [TJ It indicates that the resonance 
is not caused by the perturbation approaximation in the 
linear response approach but is attributed to the removal 
of the degeneracy of two crossing Landau levels by the 
external field. 



IV. NUMERICAL SIMULATION 

After establishing a physical picture for the electric res- 
onance of spin polarization, we come to study the effect 
of impurities. The effect can be described by introducing 
a finite lifetime t phenomenally in the Kubo formula. To 
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FIG. 3: (color online) (a) The spin susceptibility Xe as 
a function of 1/B near the resonant point for different 
electric fields at T = lCT 6 ^. (b) The electric field E 
dependence of the peak value of Xe at a low 
temperature T — 10~ 6 K. (c) The temperature T 
dependence of the peak value for an electric field 
E = l(r 6 N/C The unit of x™ is ty^l N/C. 



go beyond the linear response theory, we do numerical 
simulation to study the impurity effect in the real space. 
For this purpose, we still take the periodic boundary con- 
dition in ^-direction but an infinite potential wall in y- 
direction: V(y) = for \y\ < L a /2,and +00 otherwise. 
The disorder potential U(x, y) is taken as the short-range 
impurities of strength it, uniformly distributed at (xj, y,) 
in the plane^: 

U(x,y) =^2ui5(x-Xi)S(y-yi), (11) 

i 

where u t £ (-u/2,u/2), x { £ (~L x /2, L x /2), y { £ 
(—Ly/2, Ly/2). In the absence of impurity potential and 
Rashba coupling, the confined Landau levels \if n ks) have 
been obtained analytically 30 , which can be regarded as a 
complete set of basis. On this basis, the kinetic energy 
and the Zeeman term, has been diagonalized. The 
elements of the Rashba coupling and the disorder poten- 
tial are (<p n >k's'\H R \ip nks ) and (<p n >k's'\ U \ip„ks), respec- 
tively. After taking into account the impurity potential 
and Rashba coupling, we apply the truncation approxi- 
mation to reduce the whole Hamiltonian into an effective 
one with a finite dimension. Furthermore the effective 
Hamiltonian will be diagonalized numerically to calcu- 
late the eigenvalues and eigenfunctions. In the calcula- 
tions we take n max + 1 Landau levels n = 0, 1, ...n max 
and each Landau level has Nk discrete values of k, 
(n = 0,1,2,- • • Nk — l) 31 . Nk can also be explained 
as the maximum number of electrons accommodated in 
each Landau level. Then the number of basis functions 
we retained in the truncation approximation is N = 2x 
("max + 1) x Nk with double degeneracy of spin. By 
diagonalizing the N x N Hamiltonian numerically, one 
get A-eigenvalues E N and TV-component wave vectors 



FIG. 4: (color online) (a) Xe versus filling factor v for 

four different disorder strength A for a fixed 
temperature T = 0.025if. The inset shows the disorder 
strength dependence of the peak value, (b) Xe versus 
the filling factor v at different temperatures T for 

disorder strength A = 10 -5 . The inset is for the 
temperature T dependence of the peak value. The 
electric field is taken as low as E y = 10 -8 N/C and the 
spin susceptibility is in unit of h/4irl% N/C. 

|*m) = 5Z a nfes \^ nks ) ' wm ch are the superposition of 
the basis \<p n ks)- When a weak electric field V = e E y y is 
applied in y-direction, the expectation value of the total 
spin can be calculated numerically. As a result, Xe can 
be obtained 

L x Ly Ey 

where 

5(Sy) m = f m (^| Sy \*%)\ E ~ f m «| Sy | ^ ) | £ = Q • 

with f m the Fermi-Dirac distribution. 

In our calculations, the model parameters are taken as 
L x /N\ = Ly/Ni = y/2nlb. In this paper we take N± = 10, 
and N2 = 6. Then the maximum number of electrons at 
each level is N k = L y /Ay — N± x N 2 = 60. The num- 
ber of Landau levels are truncated to rt max = 5, hence 
N = 720. In this truncation approximation, we are con- 
cerned only with the low energy physics. The magnetic 
field is chosen as B c — 9.8T such that -Eo,+i = for 
the lowest and first excited level in the bulk region, with 
other parameters the same as in Fig. Q] 140 impurities of 
relative strength |A»| = \ui/(2irl% < A are randomly 
distributed over the sample. For each configuration of 
impurities, both the strength and position are generated 
randomly. The states arc filled from lower to higher en- 
ergy, and correspondingly, Xe can be calculated by the 
formula in Eq. (fT2]) for each configuration. 

After averaging over 10 4 different impurity configura- 
tions, we plot the average value of Xe versus the fill- 
ing factor v in Fig. 0] (a) for different strength of dis- 
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FIG. 5: (color online) The electric field E dependence 
of Xe an d the spin density S y is denoted respectively 
by the dots and diamonds in the presence of disorder. 

The temperature is set to zero. The electric spin 
susceptibility Xe i s m un it 01 fi/Anl^ N/C and the spin 
density S y is in unit of H/AttI?. 

order A. The temperature is set to T = 0.025K hence 
the ratio kT/fku c = 2 x 1CP 4 < A and the electric field 
E y = 10- & N/C such that e = eE y l b /hu c ~ 10~ 8 < A. 
When the electrons are filled upto v — 0.5, x'e displays 
a resonant peak. In contrast Xe ' IS finite in the non- 
degenerate region with the filling number v = 2 ~ 4 and 
tends to be suppressed when A = 1.6 x 10~ 3 . The rel- 
ative error due to impurity fluctuation is estimated to 
be about 0.01% and upto 5% at the resonant point due 
to the fact that the susceptibility is very sensitive to the 
energy gap opened by impurities. The peak height de- 
creases with the disorder strength A and the peak value 
versus A is demonstrated in the inset. When we extrapo- 
late to the limit A — * oo, the peak value tends to be sup- 
pressed completely. The dependence of the resonant Xe 
on the disorder strength is similar to that on the electric 
field, thus the impurity scattering opens a gap between 
the degenerate levels just like the electric field. Inversely, 
we fix the disorder strength at A = 10~ 5 and plot the 
electric spin susceptibility x'e m Fig- HI (b) at different 
temperatures with kT/huj c > A. The peak height of Xe 



decays with temperatures. The scaling behavior of the 
peak value versus T is shown in the inset. In the pres- 
ence of disorder, the dependence of the resonant peak 
of x'e on the electric field is also simulated and plotted 
in Fig. [5] denoted by the black and red dot respectively 
for the disorder strength A = 10 -5 and 10~ 4 . At low 
temperatures, when e = eE y l b /hw c <C A, Xe ' IS indepen- 
dent of the electric field E but diverges asl/E when it is 
comparable with or greater than the disorder strength A. 
The spin density S y (E) = S y (0) + x% E@ increases lin- 
early with a weaker electric field, but saturates at higher 
electric fields as plotted in Fig. [5] by the black and red 
diamond for the disorder strength A = 10~ 5 and 10~ 4 , 
respectively. 

V. SUMMARY AND DISCUSSION 

At last the occurrence of this resonance is not limited 
only in the Rashba system. The electric spin susceptibil- 
ity depends on the symmetry of spin-orbit coupling ex- 
plicitly. The present work can be generalized to a system 
with the Dresselhaus coupling, He> = I3{p x <j x — P y & y ). 
Because the Rashba coupling can be transformed to 
the Dresselhaus coupling under the transformation of 
&x — * &y, Cy — * &x and a z — * — cr^^, we conclude that it 
is x'e instead of Xe wn ich would become divergent at 
the resonant point. 

In conclusion, a tiny electric field may generate a fi- 
nite spin polarization in a disordered Rashba system in 
the presence of a magnetic field. As a result the elec- 
tric spin susceptibility exhibits a resonant peak when the 
Fermi surface goes through the crossing point of two Lan- 
dau levels. Numerical results demonstrate that the result 
goes beyond the linear response thoery. This provides a 
novel mechanism to control spin polarization efficiently 
by an electric field in semiconductors. As the spin po- 
larization can be measured very accurately it is believed 
that the effect can be verified in the samples of 2DEG, 
such as the hetero-junction of InGaAs/InAlAs. 
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